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Summary. We study the invariant measure or the stationary density of a coupled 
discrete dynamical system as a function of the coupling parameter e (0 < e < 1/4). 
The dynamical system considered is chaotic and unsynchronized for this range of 
parameter values. We find that the stationary density, restricted on the synchroniza- 
tion manifold, is a fractal function. We find the lower bound on the fractal dimension 
of the graph of this function and show that it changes continuously with the coupling 
parameter. 



1 Introduction 

Two or more coupled rhythms can under certain conditions synchronize, that 
is a definite relationship can develop between these rhythms. This phenomena 
of synchronization [1, 2] in coupled dynamical systems has acquired immense 
importance in recent years since it appears in natural phenomena as well as 
in engineering applications. What is even more interesting is the fact that 
even chaotic oscillations can synchronize. This observation is being utilized 
in secure communication [3]. Studying synchronization is also important for 
neural information processing [4, 5, 6]. 

These developments have lead to the investigation of coupled dynami- 
cal systems. The dynamical systems considered can either be continuous or 
discrete in time. Different types of couplings have been considered, like con- 
tinuously coupled or pulsed coupled. And many coupling topologies arise in 
nature as well as in human applications, like random, all-to-all, scale-free, 
small-world, nearest neighbour lattice, etc. Also the coupling strengths can 
vary from clement to element. We plan to study a system of two coupled 
maps which forms a basic building block of all these systems and allows us 
to separate the complexity due to coupling topologies from that coming from 
the chaotic nature of the dynamical system considered. As a next step, one 
can then consider various coupling topologies. We have demonstrated in [7] 
that the result of such a system of two coupled maps can be used to derive 
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the result for a globally coupled network of N maps. Here we are concerned 
with the phenomena of complete synchronization, that is, the dynamics of two 
systems becomes completely identical after the coupling parameter crosses a 
certain critical value. 

The m,odel 

We consider the following coupled map system 

X„+i = AF{Xn) := (1) 

where X = (x, y)'^ is a 2-dim column vector, A is a 2 x 2 coupling matrix and 
-F is a map from {2 = [0, 1] x [0, 1] onto itself. In the present paper we take F 
to be the extension of the tent map ft : [0, 1] [0, 1], 

(2x 0<x<l/2 
^*^^) = |2-2a; l/2<x<l ' 

to two variables and we choose 

where < e < 1 is the coupling strength. This type of coupling has been 
called contractive or dissipative. See [1] for the physical motivation behind 
considering such a system. Furthermore, the row sums of A are equal to one 
which guarantees the existence of a synchronized solution. 
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Fig. 1. Asymptotic distribution of 2000 points starting from an uniform distribution 
over the whole phase space fi. e = 0.2. 

Before proceeding we give a mathematical definition of synchronization we 
are interested in. 

Definition 1. A discrete dynamical system S : Q ^ Q is said to completely 
synchronize if as n tends to infinity |<S'"(a;) — <S'"(j/)| tends to zero for all 
{x,y) e ^. 
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Using the linear stability analysis it can be shown [8] that this coupled map 
system synchronizes when 1/4 < e < 3/4. The same result was also obtained 
by studying the evolution of the support of the invariant measure [7] which is 
a global result as opposed to the linear stability analysis which is carried out 
near the synchronized solution. This was done by showing that if e < 1/4 we 
obtain a quadrilateral of nonzero area as the support of the invariant measure. 
We depict this area in Fig 1 along with a distribution of points obtained from 
uniform distribution of initial conditions. And when e crosses the value 1/4 
this quadrilateral shrinks to the line x = y. In this sense the synchronization 
transition is a discontinuous transition. But this figure is misleading as it does 








Fig. 2. The stationary density: a histogram of the distribution of points of Fig. 1. 
Here the number of point is 10^ and the mesh size is 100x100. The horizontal plane 
depicts f2 and the vertical axis gives the number of points lying in the mesh square. 

not tell us anything about the density of points. As we show in Fig. 2, if we 
plot the histogram then we see a quite irregular structure. In Fig. 3 we plot 
the section of this invariant density along the line x = y. We see clearly that 
it is an irregular function. Studying this invariant density is the object of this 
paper. In particular we show that this density is indeed a fractal function on 
the synchronization manifold, i.e., its section along line x = y, and the fractal 
dimension of the graph of this function depends on the coupling parameter. 

The paper is organized as follows. In section 2 we give a brief introduction 
to invariant measures recalling some definitions and results required for com- 
pleteness. We also outline a method to find the stationary density. We then 
move on to our main result in the section 3 after introducing basic concepts 
needed to characterize the irregularity of fractal functions. Section 4 concludes 
by pointing out some future directions. 

2 Invariant measure 

Invariant measures or the stationary densities [9] provide a useful way to study 
the asymptotic behavior of dynamical systems. One starts with a distribution 
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Fig. 3. The section of tiio stationary density along the line x — y. In this graph the 
total number of points used is 10*^ and the size of the mesh used is 1000x1000. 



of initial conditions and studies its evolution as time goes to infinity. It is 
interesting to note that even if the dynamical system is chaotic a well behaved 
limit can exist which can then be used to study various average properties. 
In this section we give a brief introduction to invariant measures and a way 
to find one. We begin with the definitions (the norms used are norms 
throughout): 

Deflnition 2. A measure /i is said to be invariant under a transformation S 

if ^{S^^{A)) — /i(j4) for any measurable subset A of Q. 

Definition 3. Let {X,A,ijl) be a measure space and the set D[X, A, ji) be de- 
fined by D{X,A,fi) = {/ e L'^{X,A,ij,) :/> anrf ||/|| = 1}. Any function 
f G D{X, A, fi) is called a density. 

Definition 4. Let {X, A, fj) he a measure space. A linear operator P : L^ ^ 

L^ satisfying 

(a) Pf >0 for f >0, f e L^; and 

(b) \\Pf\\ = \\f\\,forf>0,feL^ 
is called a Markov operator. 

A Markov operator satisfies a contractive property, viz., ||P/|| < ||/||- And 
this property implies the stability property of iterates of Markov operators, 
viz., — Pg\\ < \\f — g\\. We shall be interested in the fixed points of 
Markov operators. 

Definition 5. Let {X,A,ii) be a measure space and P be a Markov operator. 
Any f £ D that satisfies Pf = f is called a stationary density of P. 

A stationary density is the Radon-Nikodym derivative of an invariant measure 
with respect to fx. 
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2.1 The Probenius-Perron operator 

There exist various methods to find invariant measures. One approach is to 
use the so called Probenius-Perron operator. This operator when applied to 
Pn{x), the density at the nth time step, yields the density at the (n + l)th 
time step. Since all the points at the (n + l)th step in some set D come 
from the points in the set S~^{D) we have the following equality defining the 
Probenius-Perron operator P: 

[ [ Pp{x)p{dx) = f f p{x)n{dx) (4) 

J Jd J Js-^{D) 

The Probenius-Perron operator is a Markov operator. 



2.2 The invariant measure of the tent map 

If in (4) we choose our discrete dynamical system to be the one dimensional 
tent map defined in (2) and D = [0, x] then the equation (4) reduces to 

Pp{x)=p{x/2)+p{l-x/2). (5) 

We are interested in the fixed point solutions; this leads to 

p{x)=p{x/2)+p{l-x/2). (6) 

A solution of this functional equation is p{x) = 1. Of course, this is a trivial 
solution. S{x) and S{x — 2/3) also solve this equation but we are not interested 
in such singular solutions since they do not span the phase space. 



3 Stationary density of the coupled tent map 

Wc now turn to coupled maps. In this section we study the stationary density 
on the synchronization manifold, i.e., the line x = y. We show that the density 
is a fractal function with its Holder exponent related to the coupling parameter 
e. 

We use the following definition of Holder continuity and its relation to the 
box dimension [10]: 

Definition 6. A function / : [0, 1] — > i? is in C^^, for < a < 1 and xo G 
[0,1], if for allxe [0,1] 

\f{x)-f{xo)\<c\x-xor. (7) 

A pointwise Holder exponent ap{xo) at xq is the supremum of the as for 
which the inequality (7) holds. Wc also use the relation between the Holder 
continuity and the box dimension of the graph a function, dimegraph/. 
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Proposition 1. If for f : [0, 1] — > R, the pointwise Holder exponent is a for 
all X e [0, 1] and c in (7) is uniform then diniBgraph/ = 2 — a. 

Now we are ready to state and prove our main result, viz, the estimate of the 
box dimension of the graph of the stationary density. 

Theorem 1. Let p{x,y) be the stationary density of the coupled dynamical 
system (1) and let p^{x) be its restriction on the line x = y. If p{x,y) is 
bounded then dimegraphp^ > d where d = 2 + ln(l — 2e)/ln2, with < e < 
1/4. 

Proof: We use the Frobenius-Perron operator defined in equation (4). We 
choose D = [0,x] X [0,y] and get 

Pp{x, y) = A I- /■ / p{x', y')dx'dy' (8) 
dxdyj Js-^D) 

Our S in equation (1) is not invertible. In fact, it has 4 preimages. Let us 
denote them by S^^, i = 1, ...,4. li X G O, since / is symmetric, we get 

4 

PpiX) = J-\X)J2piSr\X)) (9) 

i=l 

where J~^{X) = \dS~^{X)/dX\. The fixed point of this operator is given by 
the following functional equation for the density. 

Pix, y) = ^^^^^^^ [p(M2 - 7y/2, -IX /2 + (3y/2) 

+p{l - (3x/2 + 7?;/2, -7,t/2 + /3y/2) 

+p{Px/2 - 7y/2, 1 + 7a;/2 - Py/2) 

+p{l - l3x/2 + 72//2, 1 + 7x72 - /3y/2)] (10) 

where 7 = e/1 — 2e and /3 = 1 + 7. Since we know that a point belonging to fl 
does not leave D, all the arguments of p on the right hand side of the above 
equation should be between and 1. This gives us four linos which bound 
an area < j3x/2 — 7^/2 < 1 and < —^x/2 + j3y/2 < 1. Lets denote this 
area by Gamma. The support of the invariant measure should be contained 

in rnfi. 

We also remark that for < e < 1/4, the discrete dynamical system S 
that we have considered is everywhere expanding and this implies that the 
stationary density exists [9]. 

The above equation can be written as, for < e < 1/4, 

P{x, y) = ^ pssiPx/2 - 7y/2, -7a;/2 + fy/2) (11) 



where 
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!)) + p[l- X. ij) + (>[.)■ , 1 - ij) + (>[!- 1 - y) 
pss{x,y) = (12) 

is the part of p that is symmetric around 1/2 in both arguments. Now if we 
substitute x = y va. equation (11) we see that the arguments on both sides 
of the equation belong to the diagonal. As a result we obtain a functional 
equation 

p''{x) = ^^^p'i{xm (13) 

where wc use a shorthand notation Pg{x) = pgg{x,x). With the change of 
variable z = x/2 and a decomposition of p^{x) as p^{x) = Ps{x) + p^{x), 
the "symmetric" and "antisymmetric" part where again p^{x) is a shorthand 
notation for p^g{x, x) + Psa{x, x) + Paa{^i ^) ' ^® arrive at 

p'i{z) = {l-2e)p'^{2z)+g{z) (14) 

where g{z) = (1 — 2e)p^{2z). Its solution can be written down as 

oo 

p^(^) = ^(l-2e)^5(2'=^). (15) 
fe=o 

This is a Wcicrstrass function and if g{z) G where j3 > — ln(l — 2e)/ln2 
then the calculation in [10] for g{z) = sin(2;) can be carried over and it can be 
shown that the pointwise Holder exponent of this function is — ln(l — 2e)/ ln2 
everywhere implying that dimegraphpg = d. And ii g{z) is not smooth enough 
then it can only increase the box dimension of p^{x), hence the result. □ 



4 Concluding Remarks 

Wc have studied the stationary density of two coupled tent maps as a func- 
tion of the coupling parameter. We find that even though the density of the 
individual tent map is smooth it becomes very irregular as soon as a small 
coupling is introduced in the sense that the pointwise Holder exponent is 
small everywhere. And the density smoothes as the coupling is increased. It 
becomes smooth that is the Holder exponent becomes one at the value of 
e where the synchronization transition takes place. We have thus elucidated 
a new aspect of synchronization in coupled dynamical systems, beyond the 
standard aspects of linear or global stability of synchronized solutions. 

It is a curious fact that the Holder exponent becomes one exactly at the 
critical value of the coupling parameter, i.e., Cc = 1/4. It would be important 
to understand if there is any underlying principle behind this observation, 
that is, one valid also for other maps with varying coupling matrices. 



8 Jost and Kolwankar 



It is interesting to note that fractal probability densities have arisen in a 
completely different scenario, namely the random walk problem with shrinking 
step lengths [11]. 

One should also characterize the stationary density away from the syn- 
chronization manifold. It is expected to have a more complex multifractal 
character [12]. The effect of different network topologies on the stationary 
density is another interesting topic. 

One of us (KMK) would like to thank the Alexander-von-Humboldt- 
Stiftung for financial support. 
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